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ABSTRACT 

O 

o 

CN ' Light curves from a growing number of Gamma Ray Bursts (GRBs) indicate 

that GRBs copiously produce radioactive Ni moving outward at fractions of the 
C/^ ■ speed of light. We calculate nuclear abundances of elements accompanying the 

' outflowing Ni under the assumption that this Ni originates from a wind blown off 



(N 



of a viscous accretion disk. We also show that GRB's likely contribute appreciably 
to the galactic inventory of ^^Ca, ^^Sc, ^^Ti, ^^Ti, ^^Cu, and may be a principal 
site for the production of ^^Zn. 
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G^ ■ Subject headings: gamma rays: bursts — nucleosynthesis — accretion disks 
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^ ' 1. Introduction 



In this letter we consider the contribution of Gamma Ray Bursts to the Galactic inven- 



O 

eg ■ tory of some Fe-group elements. Our study is motivated by mounting evidence that a fair 
fraction of GRB's are associated with the production of a sizable amount of ^^Ni moving 



out with near-relativistic velocities (Price et al. 2003; Stanek et al. 2003; Hjorth et al. 2003; 
Patat et al. 2001; Iwamoto et al. 1998; Woosley & Heger 2003). The ^^Fe to which the 
radioactive Ni decays is not important for the present-day inventory of that element. As we 
discuss, however, other isotopes synthesized in the unique outflow producing ^^Ni likely are 
important for Galactic chemical evolution. 

Though there are other possibilities, we will assume that GRBs are produced by a 
viscous black-hole accretion disk formed after the collapse of a rotating massive star (Woosley 
1993; MacFadyen & Woosley 1999). Within the context of this collapsar model there are 
two possible origins for the observed Ni. As in "successful" SNe, Ni may be synthesized 
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explosively as a strong shock traverses the stellar mantle and explodes the star. Maeda & 
Nomoto (2003) have discussed nucleosynthctic consequences of this picture. Those authors 
show that a very energetic shock driven by bi-polar jets synthesizes a peculiar abundance 
pattern that may be responsible for anomalies observed in extremely metal poor stars. 

The second possibihty is that observed Ni comes from a vigorous wind blown from the 
accretion disk (MacFadyen 2003; MacFadyen & Woosley 1999). We discuss implications 
of this scenario, which is qualitatively different from shock-driven nucleosynthesis. Rather 
than develop a theory describing the disk wind we simply begin with the assumption that 
a wind is responsible for the observed Ni. As we show, this is sufficient to allow interesting 
statements about the nucleosynthesis. However, it should be kept in mind that the origin of 
the Ni and the nature of the GRB central engine remain uncertain. 

2. Nuclei Accompanying Ni in a Disk Wind 

Nucleosynthesis in the disk-wind is sensitive to the dynamic timescale rdyn characterizing 
the expansion of the wind, the entropy per baryon s in the wind, and the neutron to proton 
ratio n/p in the wind. To a large degree these parameters are constrained by observation. 
These constraints follow in part from the observations that Ni is ejected relativistically. 
Also, the efficiency for Ni production cannot be too low. This is because the disk can eject 
at most one or two solar masses of fast and initially hot material. In order to account for 
the ~ 1 /2Mq of Ni inferred from the light curves of SN1998bw and SN2003dh the mass 
fraction of Ni in the wind must satisfy X(^^Ni) > 1/4 — 1/2 if only one or two solar masses 
of material is ejected from the disk. 

The presence of any ^^Ni implies that Fg > 0.485 (Hartmann, Woosley &: El Eid 1985). 
A high efficiency for Ni production (X(^^Ni) > 0.25) implies Ye > 1/2, with the exact 
value depending on the other outflow parameters. Detailed estimates of the electron fraction 
cannot be made without reference to a speciflc disk/outflow model. This is because the 
temperature and density in the disk are high enough that weak interactions determine the 
composition of the disk (Pruet, Woosley, & Hoffman 2003; Beloborodov 2003; Surman & 
McLaughlin 2003) and the composition of the wind flowing off the disk. Pruet, Thompson 
& Hoffman (2003) employed a steady-state wind picture to model the disk outflow and found 
asymptotic electron fractions anywhere in the range 0.50 — 0.55, with the value depending 
sensitively on the accretion rate and viscosity of the disk. Surman & McLaughlin (2003) 
studied the influence of charged current neutrino capture and found that neutrino capture can 
increase the electron fraction by a considerable fraction if the neutrino luminosity of the disk 
is large. Charged current neutrino capture is expected to increase the asymptotic electron 
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fr action in outflows from a "canonical" disk with viscosity a — 0.1, mass accretion rate 
M = O.IMqScc"^ and Kerr parameter a = 0.95 by SY^. f« 0.05. For other trajectories/disk 
models the effect may be larger. The largest consistent with efficient Ni production is 
Yg PS 0.60, with the exact value again depending on other outflow parameters. 

Uncertainty in the electron fraction is somewhat of a hindrance to the present study. 

However, as we show, overall characteristics of the nucleosynthesis arc relatively insensitive 
to Yg as long as 0.505 < Fg- The special case of 0.50 < Y^. < 0.505 will be mentioned, 
but not emphasized because it seems unlikely that the electron fraction in the bulk of the 
outflow should be within 1% of the minimum value for efficient Ni production. By contrast, 
in explosive burning Ye is set not by weak interactions, but by the initial nearly isospin 
symmetric composition of the burning shell. 

To make inferences about the dynamic timescale and entropy in the wind we make the 
assumption that the wind is coasting at the low temperatures important for nucleosynthesis. 
The assumption that the wind is not accelerating at Tg = T/IO^K < 5 seems plausible since 
an accelerating wind generally expands too quickly for efficient Ni synthesis. Also there is 
little enthalpy left for driving acceleration at such low temperatures. The great kinetic energy 
of the wind and small radii at which the wind achieves low temperatures argue that outflows 
from the disk arc not decelerated significantly before nuclear recombination. However, the 
assumption of a coasting wind will have to be tested against numerical simulations that 
include expulsion of the stellar mantle in a consistent way. 

Mass conservation determines the evolution of coasting winds through M — A^r'^pVf, 
with p the density, vj the asymptotic velocity, and M an effective spherical mass outflow 
rate. For these winds the dynamic timescale is conveniently expressed in terms of the entropy 
and 



Here M_i = M/O.IMqScc"^ and Vq.i = v/O.lc. To make inferences about (3 and s, note 
that efficient Ni production implies that /3 > 3 if s = 50 and /3 > 0.1 if s = 30 (see Fig. 3 in 
Pruet, Thompson & Hoffman (2003)). Also, if the Ni outflow is to be relativistic {vq,i > 1), 
P cannot be larger than about 20. 

With estimates for Y^ and the expected range of P and s, we can calculate nucleosynthesis 
in the disk wind. In Fig. (1) we show the results of nucleosynthesis calculations for different 
assumptions about the wind parameters. These assumptions are thought to approximately 
bracket the expected range of conditions in those accretion disk outflows that efficiently 
synthesize ^^Ni. The unnormalized overproduction factor for nucleus j is defined here as 
Xj/Xqj, where Xj is the mass fraction of the nucleus j in the wind and Xqj is the mass 
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fraction of the nucleus in the sun. Though the finer details of nucleosynthesis depend on the 
wind parameters, there are some interesting features of the nucleosynthesis that do not. In 
particular, ^^Ca, ^""Sc, '^^Ti, *^^Ti, ^^Cu, and ^^Zn are all produced with large overproduction 
factors independent of the wind parameters. The abundance pattern is quite different from 
the abundance yields obtained in an explosive burning scenario (see Table 2). This is not 
surprising, since relativistic winds and explosive shock-burning are very different events. 

The influence of on the nuclear abundances is addressed in Table 1, where production 
factors for ^"''Sc and ^"^Zn are shown for a wind described by /3 = 4 and s = 30. Here the 
production factor for nucleus j is defined as 



where Mj is the total mass of the nucleus j ejected in the wind, M"'*"'^ f« 1 — 2Mq is 
the total mass coming from the disk in the form of a wind, and M^^ ~ 10 — 20Mq is the 
total mass ejected in the supernova explosion. Changing from 0.50 to 0.505 results in 
dramatic changes in the nucleosynthesis. For 0.505 < < 0.55 the production factors for 
^^Sc and ®'^Zn are within a factor of two or three of the production factors for these nuclei 
at Ye = 0.51. For larger Y^ production of ^^Sc approximately asymptotes to a maximimum 
value while production of ®^Zn decreases. As mentioned above, nuclear yields are not terribly 
sensitive to Yg unless weak interactions somehow conspire to set the electron fraction within 
1% of Ye — 0.50. As can be seen, an electron fraction very close to 1/2 might have interesting 
implications for ^^Zn. 

For definiteness we focus in the following discussion on a wind with s = 30, (3 = 4 and 
Ye = 0.51. These parameters are in-between the more extreme cases shown in Fig. (1). In 
Table 2 we show the production factors for nuclei synthesized in this wind. Also shown in 
Table II are production factors for the SN explosion of a 20Mq star as calculated by Woosley 
& Weaver (1995). Note that the production factors for the reference supernova are of order 
10, which is the typical value required to explain the presence of nuclei attributed to type II 
SNe (Mathews, Bazan, & Cowan 1992). By contrast, the production factors for the coUapsar 
wind are 100 or greater in some cases. To connect this with galactic chemical evolution, note 
that the "SN-equivalent production factor" for coUapsars is 




(2) 
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where /c is the fraction of core collapse SNe that become coUapsars. 



The fraction of SNe that become coUapsars is not well known. Observationally, the 
GRB rate is ~ 1/100 the SN rate (e.g. FraU et al. (2001)). This implies the rough lower 
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limit fc > 0.01. The true fraction is likely larger because special conditions must be met in 
order for a coUapsar to be observed as a GRB. For example, the hydrogen envelope of the 
progenitor star must have been blown off in order for a jet to make its way out of the star. 
Heger et al. (2003) estimate that fc could be as large as 0.1 depending on the metallicity 
and initial mass function for stellar formation. 

Nuclei with O'^'^^^'" ~ 5 — 10 will have significant contributions from coUapsars. In 
particular, if > 1/50, half or more of the galactic abundance of ^^Sc will come from 
collapsars. If > 0.05 - 0.1 then the abundances of ^^Ca, ^^Sc, ^^Ti, ^^Ti, siy, ^oCr, ^^Cn, 
and ®^Zn, can all have significant contributions from collapsars. A detailed discussion of 
implications of nucleosynthesis in collapsars is beyond the scope of this Letter. Here we 
discuss some interesting implications as can be inferred from the work of Woosley & Weaver 
(1995) and Timmes, Woosley & Weaver (1995). 

^^Ca: In SNe this isotope is made during explosive oxygen burning in sufficient quantities 
to explain the observed solar abundance. ^°Ca is about 100 times more abundant than ^^Ca 
in SN ejecta. In the present calculations, however, ^^Ca and ^°Ca are ejected with similar 
abundances. Future observation of a large abundance of ^^Ca relative to ^'^Ca in a metal-poor 
star might indicate enrichment by a coUapsar. 

^^Sc: This is the only stable Sc isotope. Chemical evolution studies show that at low 
metallicities Sc is under produced by about a factor of 1.5 relative to observation. Our 
calculations indicate that collapsars can synthesize sufficient Sc to explain the discrepancy 
if the electron fraction in the disk wind is larger than about 0.505. 

49rp-. rpj^-g jgQ^Qpg jg under produced by about a factor of two in chemical evolution 
studies. Collapsars can explain the discrepancy if fc ~ 1/20. 

^''Zn: The origin of this isotope is a mystery since it is under produced by about a 
factor of 5 in chemical evolution studies. One possible site for the origin of this isotope 
is the modest entropy early-time neutrino driven wind occurring after core bounce in SNe 
(Woosley & Hoffman 1992). Collapsars may be another if the electron fraction in the disk 
wind is less than about 0.55. 

In addition to implications for chemical evolution and abundances in metal-poor stars, 
the scenario we discuss also has implications for the recently observed X-ray emission lines 
that have been interpreted as iron e.g. (Piro et al. 1999; AntoneUi et al. 2000; Piro et al. 
2000), for a review see (Boettcher 2002). It has been suggested that these hues may instead 
be nickel, for example, either nickel produced in the disk that lines the falls of the hole left 
over from the jet, or small amoimts of nickel which have been ejected at rapid velocities 
behind the jet (McLaughlin, Wijers, Brown, & Bethe 2002; McLaughhn & Wijers 2002). As 
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we have discussed, if the conditions are such that ^^Ni produced in disk winds is to account 
for the observed hght curves, then it is by far the most abundantly produced element in the 
wind, with the second most abundant contributing at the 5% level or less, either ^°Zn or 
^^Zn. Therefore, if this is the mechanism for producing the observed X-ray line features, one 
expects initially a strong line or lines from the nickel, with much weaker lines from other 
elements such as zinc. Explosive burning on the other hand, produces considerable nickel, 
but also with large mass fractions of light elements such as ^^Mg. 

3. Summciry 

Observations of "''^Ni associated with GRBs hint at the occurrence of a unique nucleosyn- 
thesis event. Massive, modest entropy, rclativistic winds may be responsible for production 
of observed Ni. If so, GRBs are important for Galactic chemical evolution and can make 
important contributions to the abundances of ^^Ca, ^^Sc, ^^Ti, ^^Ti, ^^Cu, and ^^Zn. 

This work follows previous works indicating that GRBs are associated with interesting 
explosive nucleosynthesis (Maeda & Nomoto 2003). As well, GRBs may make observable 
amounts of deuterium (Pruet, Guiles, & Fuller 2002; Lemoine 2002; Beloborodov 2003), and 
may be responsible for the synthesis of some r-process elements (Pruet, Woosley, & Hoffman 
2003; Inoue et al. 2003). Gamma Ray Bursts are emerging as important sources of diverse 
nuclei. 
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Fig. 1. — Overproduction factors calculated for different accretion disk winds. Solid lines 
connect isotopes of a given element. The most abundant isotope in the Sun for a given 
element is plotted as an asterisk. A diamond around a data point indicates the production 
of that isotope as a radioactive progenitor. Here X(Ni) is the mass fraction of ^^Ni in the 
wind. Results are shown for low entropy (s/Zcb = 20) winds that inefficiently {P = 0.25) and 
efficently (/3 = 16) synthesize ^^Ni and high entropy (s = 50) winds that inefficiently (/3 = 1) 
and efficiently (/3 = 16) synthesize ^^Ni. For all cases Ye = 0.51. Note that ^^Sc, ^''Zn and 
several other Fe-group elements have large overproduction factors in all the wind models. 
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Table 1. Sensitivity of the production of ^^Zn and '^''Sc to the electron fraction in the wind' 



0.50 


3 


240 


0.505 


270 


125 


0.51 


590 


124 


0.55 


650 


33 


0.60 


590 


15 



'^Calculated assuming a total 
wind mass of IMq and a total 
stellar ejecta mass of 2OM0. 
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Table 2. Production factors for a canonical wind {P — 4:, — 0.51, s — 30) ^ 



Nucleus References O's^ 



42Ca 


260 


15 


^^Sc 


590 


5 


46Tj 


211 


11 


49Ti 


161 


7 


Sly 


61 


7 




54 


12 


53Cr 


6 


10 


55 Mn 


5 


6 


56pe 


23 


5 


59Co 


73 


7 


60Ni 


51 


5 


63Cu 


135 


26 


e^Zn 


124 


3 



^Calculated assuming a total wind 
mass of IMq and a total stellar ejecta 
mass of 20Mq. 

^Only isotopes with a production 
factor larger than 5 are shown. 

''For model S25A of Woosley & 
Weaver (1995) 



